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Abstract 

Through an integrated macroscale/mesoscale computational model, we investigate the developing shape and grain mor¬ 
phology during the melting and solidification of a weld. In addition to macroscale surface tension driven fluid flow and 
heat transfer, we predict the solidification progression using a mesoscale model accounting for realistic solidification 
kinetics, rather than quasi-equilibrium thermodynamics. The tight coupling between the macroscale and the mesoscale 
distinguishes our results from previously published studies. 

The inclusion of Marangoni driven fluid flow and heat transfer, both during heating and cooling, was found to be 
crucial for accurately predicting both weld pool shape and grain morphology. However, if only the shape of the weld 
pool is of interest, a thermodynamic quasi-equilibrium solidification model, neglecting solidification kinetics, was found 
to suffice when including fluid flow and heat transfer. 

We demonstrate that the addition of a sufficient concentration of approximately 1 pm diameter TiN grain refining 
particles effectively triggers a favorable transition from columnar dendritic to equiaxed grains, as it allows for the latter 
to heterogeneously nucleate in the undercooled melt ahead of the columnar dendritic front. This transition from columnar 
to equiaxed growth is achievable for widely differing weld conditions, and its precise nature is relatively insensitive to 
the concentration of particles and to inaccurately known model parameters. 
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1. Introduction 

Welding is a ubiquitous industrial process of great eco¬ 
nomic and technological importance [1]. Welding processes 
involve complex physical phenomena spanning multiple 
length and time scales [2]. In particular, fusion weld¬ 
ing processes involve phase changes, heat transfer by con¬ 
duction, convection and radiation, as well as surprisingly 
strong fluid flow driven by Marangoni effects (gradients in 
surface tension), all of which are tightly coupled to one 
another. 

Key mechanical properties of alloy welds are related 
both to (i) their shape and to (ii) their grain structure, and 
thus it is highly desirable to exert as much control as pos¬ 
sible on both during weld formation. Whereas the shape 
is mainly determined by macroscopic phenomena such as 
heat transfer and Marangoni driven fluid flow during the 
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melting phase, and can be largely influenced by modify¬ 
ing the fluid flow through the addition of surface active 
species [3, 4], the grain structure is mainly determined by 
mesoscopic (grain scale) crystallization phenomena during 
solidification of the melt, and can be controlled by the ad¬ 
dition of non-melting grain refining particles [5, 6, 7, 8, 9]. 

Typically, the grain morphology of a post-solidification 
weld consists of columnar dendrites, characterized by elon¬ 
gated, tree-like columns of solid which have grown into 
the melt. Equiaxed solidification, where solidification nu- 
cleation occurs within the melt away from the columnar 
front, is not common, as the thermal gradient ahead of 
the solidification front is too large to allow for the nec¬ 
essary undercooling for equiaxed growth to occur [6, 10]. 
However, the transition of the common columnar solidi¬ 
fication mode to equiaxed solidification is of practical in¬ 
terest. Equiaxed material is less prone to the unwanted 
hot-cracking (tears appearing in the fusion zone near the 
end of the solidification process [11]) and also impedes the 
undesired segregation of alloying elements to the central 
plane of the weld ]5, 7, 12]. 

Since quantitative experimental research of macroscopic 
shape evolution and mesoscopic structure evolution in 
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Nomenclature 



C 

Dendrite kinetics coefficient 

u 

Fluid velocity 

Co 

Alloy composition 

V 

Dendrite tip velocity 

Cp 

Heat capacity 

Cequi 

Volume of equiaxed dendrites 

D 

Dt 

Material derivative 

K. 

Volume of existing equiaxed dendrites 

drain 

Minimum diameter of nucleating particles 

Cnew 

Volume of newly nucleated equiaxed dendrites 

Di 

Diffusion coefficient of solute in the liquid 

Greek 

symbols 

m 

Cumulative distribution function for equiaxed 
nucleation 

a 

Partition coefficient 



e 

Emmisivity 

Cdamp 

Momentum sink term due to solidification 

V 

Laser absorptivity 

fs 

Volume fraction of solid in mush 

r 

Gibbs Thomson coefficient 

9 

Volume fraction of solid 

7 

Surface tension 

h 

Convective heat transfer coefficient 

Solid-liquid interfacial tension 

Isl 

kq 

Gaussian distribution coefficient 


Latent heat 

A 

Thermal conductivity 

L 

P 

Slope of the liquidus line 

Dynamic viscosity 

m 



M1/2 

Median inoculant diameter 

N 

Count of nucleated particles 


P 

Laser power 

4^col 

Volume fraction of columnar dendrites 

P 

Pressure 

i^equi 

Volume fraction of equiaxed dendrites after tak¬ 
ing into account grain impingement 

9 

Heat flux 

4>s 

Volume fraction of mush 

R 

Volume averaged dendrite envelop radius 

p 

Density 

Tn 

1/e^ radius for Gaussian distribution 


Inoculant particle number density 

^latent 

Pseeds 

Latent heat source term 

Stefan Boltzman constant 

CTb 

<1 

Fusion entropy 

CTd 

Standard deviation for the log-normal distribu¬ 

T 

Temperature 


tion of inoculant diameters 

at. 

Local undercooling 

Subscripts 

Too 

Ambient temperature 

n 

Normal direction 

Ts, Ti 

Solidus and liquidus temperature 

t 

Tangential direction 


welds is tremendously difficult to conduct [13], there is 
a need for efficient, realistic numerical models that can 
predict both [14]. 

The majority of previously published numerical stud¬ 
ies (e.g. [4, 15, 16, 17, 18[) on macroscopic phenomena 
during welding have focused on the evolution of the weld 
pool shape and temperature up to the end of the melt¬ 


ing stage, and thus neglect further changes in the weld 
pool shape during re-solidification after the heat source 
has been removed. However, it has been shown that the 
shape of the weld pool can still change significantly during 
the subsequent re-solidification [19, 20]. A proper study of 
macroscopic weld formation should therefore include the 
solidification stage. 
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Figure 1: Schematic representation of the studied laser welding. The 
domain is assumed to be axisymmetric. 

Other previous studies have focused on predicting the 
post-solidification mesoscopic grain structure of the weld, 
which is entirely determined during the solidification stage. 
In these studies the microstructure was either studied by 
a-posteriori analysis, neglecting undercooling of the melt 
and the possibility of heterogeneous nucleation sites (e.g. 
Zhang et al. [18]), or the solidification progression was 
modeled using mesoscopic models, while neglecting the in¬ 
fluence of macroscopic phenomena such as fluid flow in the 
molten metal on the weld pool shape [5, 21, 22, 23, 24, 25]. 

However, the macroscale and the mesoscale cannot be 
separated, as macroscale phenomena such as fluid flow and 
heat transfer determine the evolution of solidification on 
the mesoscale as well. Therefore, a proper prediction of 
weld properties requires a combination of both types of 
modeling. In this paper we present such an integrated 
macroscale/mesoscale model and we show that it can be 
used to simultaneously predict the macroscopic shape and 
the mesoscopic grain structure of a solidified conduction¬ 
mode laser weld. On the macroscale, we compute the 
heat transfer and (thermocapillary driven) fluid flow in 
the weld. On the mesoscale, the solidification evolution 
is determined from actual interface kinetics, rather than 
interface equilibrium assumptions. With this integrated 
model, we investigate the role of fluid flow during solid¬ 
ification, and the possible alteration of the solidification 
mesostructure in a steel alloy laser weld using grain refin¬ 
ing particles. 

2. Mathematical formulation 

2.1. Governing equations 

A schematic of a typical stationary weld is shown in fig¬ 
ure 1, where a non-moving slab of metal is targeted by a 
fixed high power laser. The laser irradiation will be ab¬ 
sorbed by the target material, leading to an increase in 
temperature and eventually a melting phase change. Heat 


will be transferred into the bulk of the welded material by 
conduction and convection. These phenomena are mathe¬ 
matically modeled with an energy transport equation with 
a source term for the latent heat of the phase change 

D(pCpT) = V • (AVT) + Siateut (1) 

Due to the non-uniform heating of the top surface, large 
temperature gradients will develop. These temperature 
gradients result in gradients in surface tension, leading to 
thermocapillary forces along the liquid-gas interface driv¬ 
ing flow in the weld pool. The momentum transport is 
described by the Navier-Stokes equations, with a momen¬ 
tum sink that models the friction in the so-called mushy 
zone, where the liquid and solid phase co-exist 

Df/=-Vp-I-V • (/rVH) - Fdamp (2) 

In the following sections, we will first describe the 
mesoscale models used to determine the volume fraction 
of solid and the developing mesoscale structure during so¬ 
lidification, followed by a discussion of the source terms in 
equations 1 and 2, and their boundary conditions. 

2.1.1. Columnar dendrite growth 

At the edge of the weld pool, columnar solidification 
starts instantaneously once the melt starts to be under¬ 
cooled. To compute the growth of columnar dendrites, 
we follow the model by Browne and Hunt ]26]. Here, the 
columnar dendrite front is depicted by a series of massless 
computational markers which can be thought of as an en¬ 
velope that connects all dendrite tips (see figure 2). The 
markers are displaced explicitly perpendicular to the local 
columnar dendrite front, using an analytically determined 
local growth velocity: 

V = C(Ar,)2 (3) 


8m(l — Q!)C'or 

The growth velocity v is dependent on the undercooling 
ATc = Ti — T and a kinetics coefficient C, which can be 
determined using analytical solutions of dendrite growth. 
Here, we use the model for C developed by Burden and 
Hunt ]27] (equation 4), which is based on the hypothesis 
that structures grow near the optimum condition, with Di, 
m, a, Cq and T the diffusion coefficient of solute in the liq¬ 
uid, the slope of the liquidus line, the partition coefficient, 
the alloy composition and the Gibbs Thomson coefficient, 
respectively. Alternative growth laws, such as the KGT 
model ]28] based on the marginal stability criterion ]29], 
may easily be incorporated. 

The markers are initialized along the liquidus isotherm 
once the heat source is extinguished and solidification be¬ 
gins. If a given marker temperature exceeds the liquidus 
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Figure 2: Schematic of the mushy zone at the mesoscale, with colum¬ 
nar dendrite tips, equiaxed dendritic grains and computational mark¬ 
ers (stars with velocity v) located at columnar dendrite tips 


F{d^in) = ^erfc(- ) (6) 

The number N of nucleation sites (particles larger than 
drain) ^ control volume AV is then given by 


N = PseedsAV (1 - F{dmin)) (7) 


or, in derivative form as 


dN 


\T/^F{dmm) 

P seeds 

at 


( 8 ) 


where Pseeds is the local number density of grain refiner 
seeds (seeds per unit volume). The volume of equiaxed 
grains Vequi within a control volume AV is given by the 
number of nucleated grains within AV, and the volume of 
each of those grains. The change with time in Vequi can 
then be written as 


temperature due to remelting, it is reset to the closest 
liquidus isotherm position [22]. Based on the location of 
markers, a volume fraction of mush (f>coi can be determined 
in a given finite control volume. 


2.1.2. Equiaxed nucleation and growth 

If non-melting grain refining particles are present within 
the melt, there is a possibility for nucleation of equiaxed 
grains within the weld pool away from the columnar so¬ 
lidification front. Nucleants are present because they have 
been purposely added in the form of grain-refining parti¬ 
cles to the weld pool, although in reality they may possibly 
also result from fragmentation of columnar dendrite arms. 
Recent work [30] has shown that dendrites bend / bow in 
the presence of fluid flow. Whether this leads to fragmen¬ 
tation has not been experimentally confirmed and there¬ 
fore our model takes into account added particles only. 

For the modeling of the nucleation and growth of 
equiaxed grains, we follow a volume averaging approach 
proposed by Mirihanage and Browne [31]. The minimum 
diameter of particles that can act as nucleation sites for a 
given undercooling ATc is given by 


d 


rmn 


47sz 

ASuATe 


(5) 


where ^si and ASu are the solid-liquid interfacial tension 
and the fusion entropy, respectively. 

If the grain refiner particles were equisized, we could 
assume instantaneous nucleation once the necessary un¬ 
dercooling given by equation 5 is reached. More realis¬ 
tically, however, we assume a log-normal distribution of 
grain refiner particle sizes, with a median pi /2 and a unit¬ 
less standard deviation The integral of the probability 
density function of all particles up to a size dmin gives the 
cumulative density function 


dVequi _ ^jdVex , dN 

dt ~ dt ^ dt 


(9) 


where dVexjdt is the growth of volume of existing grains 
and Vnuc is the volume of newly nucleated grains. The 
growth term can be computed as 


dVex 4^3 3 752 

— r — = -tt — R = inR V 

dt 3 dt 


( 10 ) 


The dendrite tip velocity v is computed using the same 
analytical growth velocity relationship used in the colum¬ 
nar solidification model (equation 3). The volume aver¬ 
aged dendrite envelope radius i? at a given time can easily 
be determined as 


«(*) = {''P# ( 11 ) 

From equations (5-11), we may now calculate the tempo¬ 
ral evolution of the envelope volume of equiaxed dendrites 
within the control volume AV around a location (r, AV), 
and from this we may calculate the local equiaxed grain 
volume fractions as 


/ -#. * equi 

iV) = 


Veqm{r,AV) 


AV{r) 


( 12 ) 


After some time the equiaxed grains have grown to an 
extent at which they will start feeling the presence of 
neighboring grains. To model this so-called grain impinge¬ 
ment, we restrict the growth at high volume fractions using 
an exponential relationship known as the Avrami equation 
[32]: 


(t>equi = 1 - exp(-<()*^„J (13) 
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2.1.3. Grain transport 

Unlike the columnar grains that remain stationary, 
equiaxed grains are able to move. Since we assume the 
grains perfectly follow the fluid flow, the grain transport 
can be computed by the solution of an advection equa¬ 
tion for the amount of nucleated grains and the equiaxed 
volume fraction: 


^ = v. 

dt 

{UN) 

(14) 

d(fyequi Yy 

dt 

• {Ufpequi') 

(15) 


2 . 1 . 4 . Columnar to equiaxed transition 

The columnar front progresses following the growth 
kinetics given by equation 3 until it meets a coherent 
equiaxed network. In this case, a columnar to equiaxed 
transition (GET) may take place, resulting in a different 
solidification structure. To predict this transition, for ev¬ 
ery time instance we first compute the growth of equiaxed 
grains independently of the presence of columnar grains. 
The volume fraction of the envelope of equiaxed grains is 
clipped in a way that the sum (pa = (pequi + 4>coi of equiaxed 
and columnar dendrites does not exceed 1. Now we com¬ 
pute the evolution of the columnar dendrite tips using a 
mechanical blocking criterion: We assume an equiaxed en¬ 
velope volume fraction (pg larger than a certain threshold 
to be coherent in a way that it mechanically blocks fur¬ 
ther advancement of the columnar dendrite tips. Equiaxed 
grains that are not bound in a coherent network are inte¬ 
grated into the columnar dendritic network without block¬ 
ing its progression. There is some variance in the published 
literature on the choice of the threshold value, ranging 
from 0.2 [33] to 1 [34]. Since we neglect solutal effects, 
the mechanical blocking is the only mechanism by which 
a columnar to equiaxed transition may take place. 


2.1.5. Latent heat release 

The effect of melting and solidification on the heat trans¬ 
fer are taken into account via the source term Siatent in 
equation 1. 


^latent 


= pL 


dtp 

dt 


(16) 


that envelope (the area shaded gray in figure 2) is again 
denoted by fg. Now, the solid fraction g within a control 
volume follows from 


9 = (psfs (18) 

Equation 17 has been implemented as a linear evolution 
of the solid fraction between Ti and Tg, in order to reduce 
the computational complexity of the model [26] (avoiding 
the need for iteration), as in this preliminary model the 
focus is on successfully coupling both FT and CFD models. 
A more physically realistic, non-linear relationship, e.g. 
the Lever Rule or the Scheil equation [35], could be used 
at the expense of computational efficiency. 

2.1.6. Coupling of momentum and heat transport 

Through the inclusion of the momentum sink term, the 
momentum equation 2 is valid for the entire domain in¬ 
cluding both liquid and solid regions. In the (semi-)solid 
regions, a distinction is made if it is composed of a colum¬ 
nar or equiaxed crystal structure. Columnar dendrites 
are stationary as they are attached to the unmolten solid, 
and thus are modeled as a porous medium, introducing 
a momentum sink following the isotropic Blake-Kozeny 
model [36] 


damp — 



(19) 


K = Ko 


{1-gr + e 


( 20 ) 


with pKo = 10®Ns and e = 10“®. 

Equiaxed grains are able to move with the liquid, and 
with increasing equiaxed volume fraction the liquid metal 
will turn into a slurry. We model this by increasing the 
viscosity in the slurry, following a correlation suggested by 
Thomas [37] for spherical solid particles in a liquid: 


l^equi — M 


1 -|- 2.5pequifs + 10.05(^e9«i/s)^ + 


0.00273 exp(16.6())eq„i/s) 


( 21 ) 


with g the volume fraction of solid material. While the 
heat source is active, g is evaluated as the equilibrium 
change of volume fraction of solid fg 

fg = ^l^,Tg<T<Ti (17) 

During the solidification stage, when the heat source is 
switched off, the solid fraction g is determined from non¬ 
equilibrium solidification kinetics as a combination of the 
growth of equiaxed and columnar dendrites. The envelope 
fraction (volume fraction of mush) within a control volume 
is now given by pg = pcoi+pequi- The solid fraction within 


This model breaks down when the equiaxed network be¬ 
comes coherent at high volume fractions. Thus, we switch 
to the porous medium model for equiaxed regions with a 
high volume fraction above the coherency threshold. 

2.2. Boundary eonditions 

We assume the weld pool to be axisymmetrical and make 
use of this by only simulating a wedge of the domain. Cir¬ 
cumferential gradients and velocities are zero on the wedge 
faces. The boundary conditions for the other faces are out¬ 
lined in the following. 


5 



2.2.1. Heat input 

At the top surface, the laser irradiation is modeled by a 
Gaussian distributed heat flux. At all surfaces, including 
the top surface, there is an outflux of heat due to nat¬ 
ural convection and radiation. This heat outflux is very 
small compared to the laser irradiation, but once the lat¬ 
ter is switched off, it is this heat outflux which together 
with conduction of excess heat into the unmolten material 
causes solidification of the weld. 


dT 

— Qlaser Eradiation Econvection 

(22) 

Eradiation 

= ake{T^-T^) 

(22a) 

Qconvection 

1 

II 

(22b) 

Qlaser 

= kq^exp{-kq^) 

q ' q 

(22c) 


The numerical values of the input parameters are given 
in table 2. 

2.2.2. Momentum 

At the top gas-liquid interface, we introduce a shear 
stress in the liquid due to surface tension gradients (i.e. 
Marangoni forces) along the interface: 

MV„C/i = (23) 

The variation of surface tension with temperature is 
computed using the thermochemical model of Sahoo et al. 
[38]. 

At all other surfaces, we set the velocity to zero. 

3. Numerical procedure 

Our solver is built on top of the open source finite vol¬ 
ume framework OpenFOAM (version 2.1.x) [39]. Postpro¬ 
cessing of the results for publication is done with the open 
source software matplotlib [40]. 

The non-linearity associated with the pressure-velocity¬ 
coupling is handled by the iterative PISO algorithm [41]. 
Once a divergence free velocity held has been computed, 
the temperature equation is solved. If a phase change 
occurs, the temperature equation will be non-linear and 
thus solved iteratively. The non-linearity due to latent 
heat during melting is dealt with using an implicit source 
term linearization technique [42]. Here, we assume instan¬ 
taneous progression of the melting front based on local 
temperature conditions. While this algorithm is appropri¬ 
ate for melting conditions, during solidification it fails to 
predict undercooling and the metallurgically relevant tran¬ 
sition from columnar to equiaxed grain growth. This de¬ 
ficiency is alleviated by using the mesoscale front-tracking 


model for solidification predictions by Browne and cowork¬ 
ers [26]. In this case we displace the computational mark¬ 
ers forming the solidification front based on the temper¬ 
ature of a previous time step, and then iteratively deter¬ 
mine the current temperature and solid fraction within the 
mushy zone (/^ in equation 17) with under-relaxation. 

Since OpenFOAM can perform 3D simulations only, 
the solution domain is an axisymmetric wedge of size 
r = 10 mm and z = 5 mm with an opening angle of 5°. The 
axisymmetric model presents a marked improvement over 
previously published 2D planar results e.g. by Duggan 
et al. [22] or Koseki et al. [24]. For the fluid flow and heat 
transfer equations, we use a mesh of 200x100 uniform cells. 
For the volume averaged approach to determine equiaxed 
nucleation and growth we use a mesh that is coarse enough 
to ensure a reasonable amount of grain refining particles 
per control volume. Here, we use a uniform mesh of 160x80 
cells, corresponding to an average of four grain refining 
particles per grid cell at Pseeds = 1000^/^mm The code 
is arranged in such a way that the mesh for the fluid flow 
and heat transfer calculations can be chosen independently 
from the mesh used for the nucleation and growth calcu¬ 
lations. The computations on the two meshes are fully 
coupled, meaning that data (temperature from macro- to 
meso-, and solid fraction from meso- to macroscale) is ex¬ 
changed in both directions at every time step via a linear 
interpolation framework [2[. This allows for an indepen¬ 
dent adaptation of the fluid flow and heat transfer mesh to 
the mesh requirements imposed by the macroscopic length 
scales. 

The time step is fixed at 25 ps. We use a 2nd order back¬ 
ward differencing time marching scheme, and a 2nd order 
TVD scheme (limitedLinear) for the divergence terms. 

4. Results and discussion 

Following the previous studies by Villafuerte et al. [5, 7], 
Park [9], Duggan et al. [22] and Koseki et al. [24, 43], we in¬ 
vestigate a laser spot weld on an Fe-15.9wt%-Cr-14.1wt%- 
Ni steel alloy. To determine the growth kinetics coefficient 
C, this ternary alloy is treated as a pseudo-binary alloy, as 
the chromium segregates preferentially to the liquid and 
the partition coefficient of the nickel is close to 1.0, and 
thus the chromium is assumed to be the dominant solute 

[44]. 

In contrast to previously published computational re¬ 
sults, we include fluid flow in our analysis, as it is known 
to have a strong effect on the heat transfer and shape evo¬ 
lution in a weld pool, both during the melting stage and 
the solidification stage. The influence of fluid flow on the 
mesoscale grain structure formation during solidification 
in welding has not been reported in the literature to date. 

4-.!. Weld evolution during the melting stage 

The fluid flow in the pool is driven by surface tension 
gradients due to temperature gradients. The surface ten¬ 
sion dependency on temperature is tightly coupled to the 
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concentration of surface active agents (surfactants) present 
in any weld pool, such as sulfur and oxygen. In the follow¬ 
ing we neglect the presence of oxygen and assume sulfur 
is the only surfactant (homogeneously) present in the sys¬ 
tem. 

A sketch of the problem investigated is shown in fig¬ 
ure 1. The relevant material properties together with their 
respective literature sources are listed in table 1. The prop¬ 
erties related to heat inputs and losses are listed in table 2. 
For simulations with grain refining particles present, un¬ 
less otherwise specified we use a TiN particle density of 
(1000/mm^)^/^, which is the 3D equivalent of a 2D density 
of 1000/mm^ as used by Koseki et al. [24]. In their simula¬ 
tions, Koseki et al. assign a homogeneous fixed size of grain 
refining particles based on a measured required undercool¬ 
ing of ATc = 1.8 K by Bramfitt [8] to achieve nucleation of 
equiaxed grains on TiN particles. Here, we assume a log¬ 
normal distribution of particle sizes (equation 6) as mea¬ 
sured by Park [9] with parameters /ii /2 = 0.636 x 10“®m 
and ad = 0.1, such that the peak of the distribution is close 
to the value determined by Bramfitt, while the spread of 
the distribution lies in-between the fixed size assumed by 
Koseki et al. and the wider log-normal distribution mea¬ 
sured by Park. 

The high power laser irradiation will lead to melting 
of the base metal in fractions of a second and thus the 
formation of a weld pool. For an alloy with a sulfur con¬ 
centration of 200 ppm the surface tension gradient dj/dT 
is positive for temperatures up to nearly 2100K (figure 3). 
Since fluids flow towards regions of higher surface ten¬ 
sion, the fluid flow will dominantly be directed towards 
the center of the pool, where the highest temperatures 
lead to the highest surface tension. Such a flow pattern 
leads to a weld pool with a width over depth ratio smaller 
than 1, characteristic for conduction mode welds (a weld 
with negligible vaporization of the melt [46]) with a high 
content of surface active species (figure 4, left column). 
At lower sulfur concentrations, the sign change of dj/dT 
is shifted towards lower temperatures (figure 3). For a 
concentration of 80 ppm, dj/dT will not remain positive 
over the entire weld pool surface, and thus the maximum 
surface tension will not be located at the center of the 
weld pool but shifted radially outwards, close to the crit¬ 
ical 1900 K isotherm. The liquid metal in the weld pool 
will flow towards this point, leading to two recirculation 
zones in the weld pool, resulting in decreased penetration 
with a characteristic bulge of the pool boundary under 
the stagnation point and a width to depth ratio larger 
than 1 (figure 4, middle column). The computed maxi¬ 
mum velocities in the weld pool are in the order of 0.3- 
0.4 ms“^, which is in good agreement with a maximum 
velocity of [{dj/dT)q/{^Cp)]°'^ = 0.3ms“^ obtained by a 
scaling analysis [47], and simulation results of 0.4 ms“^ for 
comparable conduction mode laser welds [4[. 

Here it is illustrative to stress that neglecting fluid flow 
entirely, as has been done in previous studies [5, 22, 23, 
24, 25], leads to a very different, hemispherical weld pool 



T [K] 

Figure 3: Change of surface tension temperature coefficient with 
temperature and impurity concentration 

shape, irrespective of the sulfur concentration (figure 4, 
right column). This highlights the need of a proper weld 
pool flow model to obtain an accurate initial condition for 
the solidification computations. 

4-2. Weld evolution during the solidification stage 

Whereas the final shape of the weld is mostly (but not 
entirely) determined during the heating phase, its grain 
morphology is fully determined during the cooling phase. 
Here we study the influence of accounting for solidifica¬ 
tion kinetics using our front tracking model versus the use 
of an equilibrium enthalpy method, and the influence of 
accounting for fluid flow, on the shape and morphology 
evolution during the cooling stage. 

For a 200 ppm sulfur weld, the evolution of the melt pool 
and temperatures therein, at two time instances after deac¬ 
tivation of the heat source, is shown in figure 5, computed 
both with the commonly used enthalpy method which does 
not take into account growth kinetics, and with our com¬ 
bined front-tracking equiaxed growth model. With the 
latter the temperatures within the weld pool show some 
difference when compared to the computed temperatures 
using an enthalpy method. Shown are the Tg = 1679 K 
solidus temperature isotherm and the Ti = 1710 K liquidus 
isotherm. Notable differences are found in the location of 
the liquidus isotherm and these differences increase dur¬ 
ing further solidification. In the enthalpy method the so¬ 
lidification front is coincident with the liquidus isotherm, 
whereas in the front-tracking method the solidification 
front progression is significantly delayed resulting in a large 
undercooled region and a much thinner mushy zone. The 
same conclusion can be drawn for the completely differ¬ 
ently shaped weld in a low (80 ppm) sulfur alloy (see fig¬ 
ure 6). The transient position of the solidus isotherm with 
the two models, on the other hand, is virtually the same 
for both the high and low sulfur case. 
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Figure 4: Development of the weld pool with isotherms in Kelvin and flow velocity vectors during the heating stage at three time instances with 200 ppm sulfur (left column), 80 ppm 
sulfur (middle column) and pure diffusion. The final image shows the state at t=0.8s when the laser is shut off and solidification begins. 
































Table 1: Material properties of the Fe-Cr-Ni alloy 


Property 

Value 

Unit 

Ref. 

Eutectic temperature T, 

1679 

K 

[22] 

Liquidus temperature T; 

1710 

K 

[22] 

Specific heat capacity Cp 

780 

Jkg-iR-i 

[22] 

Density p 

7250 

kgm“® 

[22] 

Thermal conductivity A 

35 

Wm-^K-i 

[22] 

Latent heat of fusion L 

1.854 • 10® 

Jkg-i 

[22] 

Dynamic viscosity p 

6.1625 • 10“® 

Pas 

[45] 

Surface tension temperature coefficient d^jdT 

1 

o 

CO 

1 

Nm-^K"^ 

[38] 

Entropy factor 

CO 

1 

o 

00 

CO 

- 

[38] 

Entropy of segregation 

-1.66 • 10® 

Jkg-i 

[38] 

Surface excess at saturation 

1.3 • 10“® 

kmol m~^ 

[38] 

Burden and Hunt growth kinetics coefficient C 

8.995- 10“® 

ms“^K-2 

[35] 

Diffusivity of Cr in the liquid Di 

1.8 • 10“® 

m^s-i 

[35] 

Gibbs-Thomson coefficient P 

3.88- 10"^ 

m K 

[35] 

Liquidus slope m 

-2.9 

K/wt% 

[35] 

Partition coefficient a 

0.86 

- 

[35] 

Solid-liquid interfacial tension 

0.299 

J m~^ 

[22] 

Fusion entropy 

7.71 ■ 10® 

Jkg-iR-i 

[22] 


Table 2: Properties related to the laser and heat losses 


Property 

Value 

Unit 

Laser power P 

1000 

W 

Distribution coefficient kq 

2 

[-] 

Beam radius Vq 

2.5 • 10"® 

m 

Incidence time 

0.8 

s 

Absorptivity g 

1 

[-1 

Emissivity e 

0.5 

[-1 

Convective heat transfer 

coefficient h 

8 

Wm"® K-i 

Ambient temperature Too 

300 

K 


Thus, simulations aimed solely at predicting the weld 
pool shape can be carried out with a basic enthalpy 
method without the need for a sophisticated solidification 
model. If the post-solidification grain structure is of inter¬ 
est however, a sophisticated solidification model, such as 
the one presented here, is necessary to compute accurate 
thermal gradients and solidification front progression, as 
enthalpy methods have no microstructure information. 

The question of the significance of fluid flow during the 
cooling phase on the evolution of the solidification process 
is addressed in figures 7 and 8, for high and low sulfur 
concentrations respectively. For both cases, starting from 
the same initial condition at the end of the melting phase, 
the solidification phase is subsequently computed with and 
without fluid flow. The inclusion of fluid flow in the so¬ 
lidification phase leads to better mixing and heat transfer, 
and thus lower and more uniform temperatures in the core 
of the weld pool. This leads to significant differences in 
the weld pool shape [19, 20], and more importantly in the 


post-solidification grain morphology, as will be shown in 
the next section. 

For the high sulfur case, the inclusion of fluid flow, which 
in this case is directed inward along the weld top surface 
and downward along the weld axis, leads to continued heat 
transport to the bottom of the weld after the heat source 
has been switched off, and thus to continued local melting. 
This causes a deeper and less wide weld as compared to 
the situation in which fluid flow is ignored. 

For the low sulfur case, the downward fluid flow and 
convective heat transport is weaker and further away from 
the axis. Yet, the inclusion of fluid flow leads to contin¬ 
ued heat transport to the bottom of the weld after the 
heat source has been switched off, and a deeper weld as 
compared to the situation in which fluid flow is ignored. 

4-.3. Influence of grain refining particles on grain morphol¬ 
ogy evolution during solidification 
Up to this point, we have investigated solidification 
within the weld pool in the absence of grain refining par¬ 
ticles, leading to columnar dendritic growth only. When 
taking into account the promotion of equiaxed nucleation 
using grain refiners, at some point the columnar front pro¬ 
gression is blocked and the solidification structure transi¬ 
tions into equiaxed grains. This is addressed in the current 
subsection. 

4-3.1. Sensitivity to numerical and model parameters 
First, we assess the sensitivity of the main simulation 
outcome of interest, viz. the location of the columnar- 
equiaxed transition line, as represented by the total volume 
of equiaxed solid in the weld, on (i) time step and mesh 
size and (ii) the precise value of the threshold for coherence 
of equiaxed grains. 
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Figure 5: Temperature evolution at two time instances after the 

heat source has been switched off at t=0.80s, for 200 ppm sulfur Figure 6: Temperature evolution at two time instances after the heat 

concentration. The location of the columnar dendritic solidification source has been switched off, 80 ppm sulfur 

front is depicted by a line with triangle symbols. 
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Figure 7: Temperature evolution at two time instances during solid¬ 
ification with and without the consideration of fluid flow, 200 ppm 
sulfur. The line with triangle symbols depicts the location of the 
columnar dendritic solidification front. 


Figure 8: Temperature evolution at two time instances during so¬ 
lidification with and without the consideration of fluid flow, 80 ppm 
sulfur. The line with triangle symbols depicts the location of the 
columnar dendritic solidification front. 
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Figure 9: Sensitivity on the coherency threshold of the computed 
transition from columnar to equiaxed solidification structure in a 
weld pool, 200 ppm sulfur 

Ad (i): A simulation with a doubled fixed time step of 
50 ps resulted in a slight increase of 4.6% of the volume of 
equiaxed solid (and a corresponding decrease in the vol¬ 
ume of columnar dendritic solid). The total volumes of 
equiaxed solid in the weld, as computed with a fixed time 
step of 25 ps on three different meshes, viz. the standard 
mesh, and two meshes that are coarsened and refined by 
a linear factor 1.25 compared to the standard mesh (see 
Roache [48] for a justification of the refinement ratio), 
show a monotonic convergence, allowing for the computa¬ 
tion of an asymptotic value by Richardson extrapolation. 
Compared to the asymptotic value, the coarse mesh result 
is 6.1% smaller, the base mesh result 2.0% smaller and the 
fine mesh result 0.48% smaller. 

Ad (ii): In a previous section, we mentioned that the 
precise value of the threshold for coherence of the equiaxed 
grains is not agreed upon in the literature (ranging from 
0.2 [33] to 1 [34]). To study the sensitivity of our results to 
this threshold value, we simulated solidification of a weld 
pool using a coherency threshold of 0.49 and 0.6, leading to 
a 15% difference in the volume of equiaxed structure in the 
weld (figure 9). The same conclusion has been drawn by 
Mirihanage and Browne [31], who found the occurrence of 
the columnar-to-equiaxed transition (GET) in a casting to 
be rather insensitive to the exact value of the coherency 
threshold. All further simulations in the current study 
were performed with the most commonly used blocking 
threshold of 0.49 as proposed by Hunt [49]. 

4-3.2. Influence of grain transport 

The presence of fluid flow in the weld pool during so¬ 
lidification leads to a significantly earlier transition to 
equiaxed growth, as can be seen in figures 10 and 11 for 
high and low sulfur cases respectively. Conversely, the lo¬ 
cation of the columnar-to-equiaxed transition was found 
to be rather insensitive to the inclusion or neglect of the 



Figure 10: Equiaxed transition in a weld pool with and without fluid 
flow, 200 ppm sulfur 



r [mm] 


Figure 11: Equiaxed transition in a weld pool with and without fluid 
flow, 80 ppm sulfur 

transport of equiaxed grains with the flow, as described by 
equations 14 and 15. 

Hence, the earlier transition to equiaxed growth when 
including fluid flow is mainly due to the lower thermal gra¬ 
dients and lower temperatures in the melt (see figure 7), 
which result in a larger liquid region with sufficient un¬ 
dercooling for the equiaxed nucleation and growth to take 
place. This conclusion agrees with Hunt’s criterion [49], 
stating that the fraction of equiaxed grains is inversely 
proportional to the magnitude of the thermal gradients 
cubed in the melt. 

4-3.3. Influence of grain refining particle number density 

The influence of the number density of grain refining 
particles on the transition from columnar to equiaxed so¬ 
lidification is illustrated in figures 12 and 13, which show 
that the transition to equiaxed solidification can be trig¬ 
gered with a grain refiner density of 500^/^ particles per 
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Figure 12: Equiaxed transition in a weld pool with two different 
densities of grain refining particles, 200 ppm sulfur 

mm^, and that the position of the transition is only slightly 
affected by an increase of the grain refiner number density 
to 1000^/^ particles per mm^. This observation matches 
the findings by Koseki et al. [24] that grain refiner den¬ 
sity mainly influences the size of the resulting equiaxed 
grains, and not the point at which they reach coherency. 
Simulating the effect of a grain refiner density lower than 
500^/^ particles per mm^ is not feasible using our cur¬ 
rent volume averaged approach, as this would require very 
large computational cells in relation to the weld pool area 
to allow for a sufficient number of grain refiner particles 
per computational cell. In figures 14 and 15, we show the 
thermal gradient and equiaxed envelope volume along a di¬ 
agonal line rz = + z^, z = r at a time instance where 

the equiaxed grains reach coherency. We can identify a 
low, fairly uniform temperature gradient within the weld 
pool which rises sharply around the liquid-solid interface 
and is insensitive to the number of grain refining particles 
present. As a result, favorable conditions for the growth of 
equiaxed grains are met at the same location for both par¬ 
ticle densities, though coherency is achieved slightly later 
in time for the lower number density. 

5. Conclusion 

In addition to the macroscopic shape of a weld, the 
mesoscopic grain morphology is crucial in determining the 
resultant properties. These two are determined by the 
interaction between macroscale fluid flow and heat trans¬ 
fer, and mesoscale grain solidification kinetics, during the 
heating and cooling phase of a welding process. Whereas 
both macroscale and mesoscale phenomena have previ¬ 
ously been studied separately through computational mod¬ 
els, an integrated model for both scales is needed to predict 
weld properties on both relevant scales. 

Such an integrated model has been presented in the cur¬ 
rent work. The model is conceptually simple and com- 



Figure 13: Equiaxed transition in a weld pool with two different 
densities of grain refining particles, 80 ppm sulfur 



Figure 14: Thermal gradient and equiaxed envelope volume fraction 
along a diagonal line z = r away from the weld pool center, 200 ppm 



Figure 15: Thermal gradient and equiaxed envelope volume fraction 
along a diagonal line z = r away from the weld pool center, 80 ppm 
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putationally inexpensive, while granting insight into the 
post-solidification shape and grain morphology of a weld. 

The model has been applied to study the melting and 
solidification of a laser spot weld on two Fe-Cr-Ni steel al¬ 
loys, viz. with a low and high sulfur content. Due to their 
very different dependency of surface tension on tempera¬ 
ture, the weld pools of these two alloys represent the two 
extremes in surface tension driven flow commonly observed 
during welding. Although demonstrated here for laser spot 
welding, the presented results are relevant for arc welding 
as well, as both types of welding are equivalent once the 
heat source is turned off. 

The predicted solidification evolution has been com¬ 
pared to that obtained with the commonly used equilib¬ 
rium enthalpy method, in which solidification takes place 
quasi-statically and its kinetics is ignored. Apart from the 
obvious fact that an equilibrium method cannot predict 
the important transitions in grain morphology, a marked 
difference between the two methods was observed in the 
thermal gradients within the weld pool during solidifica¬ 
tion. This in turn has a strong impact on the developing 
grain structure. However, if only the weld pool shape, and 
not its grain morphology, is of interest, the equilibrium en¬ 
thalpy method is found to be sufficient to obtain accurate 
results for the partial penetration cases examined here. 

Neglecting fluid flow during the cooling phase, as done 
in previously reported studies, was found to have a signifi¬ 
cant influence on the predicted weld pool shape and grain 
morphology. Thus any accurate simulation of weld pools 
should include fluid flow. 

The use of grain refining particles was found to be an 
effective means to favorably alter the grain morphology by 
initiating a transition from columnar to equiaxed growth. 
This transition was found to be not very sensitive to the 
most uncertain parameter in the model, viz. the coherency 
threshold value, nor to the density of grain refining parti¬ 
cles above a certain critical value. 
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